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The method of optimized perturbation theory (OPT) is used to study the phase diagram of the 
massless Gross-Neveu model in 2 + 1 dimensions. In the temperature and chemical potential plane, 
our results give strong support to the existence of a tricritical point and line of first order phase 
transition, previously only suspected to exist from extensive lattice Monte Carlo simulations. In 
addition of presenting these results we discuss how the OPT can be implemented in conjunction 
. with the Landau expansion in order to determine all the relevant critical quantities. 
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I. INTRODUCTION 



' The use of exactly solvable models and methods leading to analytical solutions has always been important in 
physics in general. Models of particular interest are those that display some physical correspondence to the physics 
^ • of quantum chromodynamics (QCD) or that may give insights to the difficult task, given its nonperturbative nature, 
^ ' of analyzing the phase structure at finite temperatures and densities. This is also the kind of situation, regarding 
, phase transitions in a hot and dense medium, that is expected to be realized in the early universe and in the recent 
experiments involving heavy-ion collisions. The use of appropriate techniques and models in the context of quantum 
, field theory that can increase our knowledge on this kind of phenomena are then timely and important. One such 
^ ■ phenomenon of topical interest is the study of the chiral symmetry restoration (CSR) in QCD. However, due to the 
intrinsic difficulty of approaching directly this study in the context of QCD itself, other more tractable models have 
■ been extensively used. Among them, the simplest fermionic model displaying many properties common to QCD 
' is the Gross-Neveu (GN) model p|. For instance, in 2-1-1 dimensions it can exhibit a discrete spontaneous chiral 
, symmetry breaking (CSB) that happens for a critical coupling at zero temperature and density, while its spectrum 
I/"} can contain excitations like baryons as well as mesons (which are composite fermionic states). Even though it is a 
, quartic interacting fermion model, thus nonrenormalizable (in 2-1-1 dimensions) in ordinary perturbation theory, it 
is renormalizable in a 1/N expansion (where N denotes the number of fermion species) and it is exactly soluble in 
the large- limit These properties make this effective model one of the simplest and most useful frameworks for 
understanding some of the important physical features displayed by QCD. 

For the reasons above, the GN model is a perfect laboratory for studying symmetry aspects in a hot and dense 
^ medium and outside the QCD realm, the model by its own and its extensions may also have applications in condensed 
matter physics in the understanding of the properties of fermion liquids and high-Tc superconductivity models 
Besides being studied analytically mostly with the large- TV expansion technique, there are a large number of works 
studying the GN model numerically from the lattice Monte-Carlo perspective 0, Q (for an introduction to lattice 
QCD, its problems in implementation and applications to effective models, like the GN one, see [6| and references 
therein) . Among the interesting thermodynamical properties of the GN model demonstrated by the large- A^ technique 
is the existence of a line of second order CSR in the temperature and chemical potential (T, /i) plane, starting at a 
value (Tc, 0) and ending, abruptly, at a first order transition point (0, /Xe) While this phase structure was verified 

by lattice simulations with finite number of fermion species, A'^, in Ref. [5| some indication was given for the existence 
of a first order transition line starting at a tricritical point (Ttcr,Mtcr) and ending at (0,/ic). However, the precise 
position of this tricritical point in the phase diagram could not be determined. Even within the 1/A^ expansion, there 
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are no results available on the phase structure of the 2+1 dimensional GN model beyond the leading order, due to 
the appearance of complicate integrals (note however that there are results up to O^l/N"^) for the critical exponents 
of the bulk transition). It then becomes an important and interesting issue to investigate whether there is indeed 
a tricritical point for the GN model at finite N in 2+1 dimensions and by which method we would be able to localize 
it in a robust form. By doing so, we expect to have a much better understanding of the chiral symmetry restoration 
problem at finite T and /i not only in this model but also to have eventually a better knowledge of a similar behavior 
that may happen in QCD or even in analog condensed matter systems. In this work we study the 2+1 dimensional 
GN model beyond the large- limit and, to obtain its phase diagram at finite T and /x, we make use of the method of 
the optimized perturbation theory (OPT) (also known as linear 6 expansion) Q. Here, we discuss how this method 
can be used to extract useful information, regarding critical quantities, from the Landau expansion of the free energy 
density. 

This paper is organized as follows. In the next section we introduce the model, presenting the OPT free energy 
density. In Sec. Ill, we discuss optimization issues related to the Landau expansion of the free energy density. Our 
results for the (T, /i) phase diagram as well as for the critical quantities are presented in Sec. IV, showing the existence 
of tricritical point for finite values of TV. Our conclusions and perspectives for future applications are discussed in 
Sec. V. 



II. THE INTERPOLATED GROSS-NEVEU MODEL 



The GN model is described by a Lagrangian density for a N component fermion field 'ipj (j = 1, . . . , N) with local 
quartic fermion interaction which is defined by [jj (in this work all our expressions are defined in Euclidean space) 



= ^i^3-i^{i>,i^,? , (1) 

where, in 2+1 dimensions, the fermion fields are treated as four-component Dirac spinors and the 7^ matrices are 
defined in the 4x4 representation. In this way, Eq. ([1]) is invariant under the discrete chiral transformations ■(/'j — * TsV'i 
and — > —ipjjz. Though it also has a continuous flavor symmetry we will only be interested in the case of discrete 
symmetry, recalling a general no-go theorem in 2+1 dimensions that forbids spontaneous symmetry breaking of a 
continuous symmetry at any flnite temperature [lol | . 

The implementation of the OPT within this model is reviewed in a previous application 11] to the 1+1 dimensional 
case. In practice, one considers the original theory, Eq. ([T]), adding a quadratic (in the fermion field) term (l — 5)'q%pjtpj 
where 77 is an arbitrary mass parameter and 5 is a bookkeeping parameter which labels the order at which the OPT 
is being performed. At the same time, A 5\. Also, as conventional, we rewrite the quartic term by introducing an 
auxiliary scalar field a. Then, Eq. ([T]) acquires its OPT form [ll[ 



Cs = 'i/'j ^tpj + Satpjipj + (1 - (5)77V'j + ^'^^ • (2) 

The OPT then describes an interpolated theory with a fermion propagator [/^ + (1 — S)ri + 6ac]~^, Yukawa vertex 
S and a propagator X/{NS). Any quantity ^^'^^ is then evaluated to an order-J*^ using these new Feynman rules. 
As shown in the many previous applications using this method jT^], it operates for any iV and is free from infrared 
divergences. Nonperturbative results for the original (ry-independent) theory are obtained by optimizing with 
respect to rj at 6 = 1. One standard and largely adopted optimization procedure is the variational criterion known as 
the Principle of Minimal Sensitivity (PMS) which fixes 1] by requiring, order by order, that d^^''^ /drj = (T3| . 

From Eq. we are interested in evaluating the effective potential for the auxiliary scalar background field, Ves{ac), 
obtained once all fermion fields and fluctuations around the scalar background field are integrated out 0- From this 
quantity all the thermodynamics for the model can be derived. In particular, CSB is signaled by a non-vanishing 
vacuum expectation value for tr, CTc = (u), which is a minimum of Vcsio'c)- 

In Fig. [T]we show all the Feynman diagrams contributing to T4ff((Tc) up to 0{6). Note that already at this first 
order in the OPT a correction that goes beyond the large- expansion is included. The first two diagrams in Fig. [1] 
are 0{1/N'^), while the third one contains an exchange type of self-energy as insertion and brings the first 0{1/N) 
correction to Vcs {(Jc). Explicitly, we have for all the terms shown in Fig. [T]the expression 
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FIG. 1: Feynman diagrams contributing to Vcff to 0{5^). The continuous lines are (the 5 independent) fermion propagators, 
the black dot is an insertion of S{ac — rf), while the dashed line represents the a propagator. 

where the EucHdean momenta are given by P = (PqiP) and Pq — {2n + 1)ttT — ifj,, with n — 0, ±1, ±2, .... In the 
above equation we have used the compact notation: 
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with all momentum space integrals done with dimensional regularization, d = 3 — e, while rcnormalization can be 
implemented in the modified minimal subtraction scheme (MS). All the contributions considered in Eq. ([3]) are 
actually finite in dimensional regularization. In Ref. [T^ we have considered the free energy up to order-(5^, at T = 
and fi — 0, where the first divergences show up and we refer the interested reader to that work for a detailed account 
on the renormalization problem of the 2+1 dimensional GN model in the OPT. Let us only emphasize here that 
though the model is not perturbatively renormalizable, to 0{6-) and within dimensional regularization, we only need 
standard (i.e. mass, wave-function, etc) renormalization counterterms to render the effective potential calculation 
finite [ij. 

III. LANDAU'S EXPANSION AND THE OPTIMIZATION PROCEDURE 

The effective potential, VcsicTc), is a natural quantity to be optimized in the spirit of the PMS procedure since it 
represents the system's (Landau) free energy density, whose knowledge allows for a complete analysis of all relevant 



thermodynamical properties. Also, as explained in Ref. by optimizing the free energy we can guarantee that not 



only the exact large-iV results are already obtained at leading order in the OPT, but also at any subsequent order, 
provided that one stays within the N ^ oo limit (at first order in the OPT this corresponds to dropping the third 
diagram in Fig. [1] and applying the optimization procedure to the first two diagrams) . This shows that the OPT 
is convergent at any order in the OPT when contrasted with the large-iV predictions. This result is valid for any 
temperature or value of chemical potential and applies as well to any number of dimensions (for other accounts on 
the convergence of the OPT in the context of quantum field theory, please see [l5| and references therein). Let us 
consider the expression for Vcf[{cTc) at 0{S) with A — > —A in which case chiral symmetry is broken at T = = /x. 
After taking the traces in Eq. ([3]), rearranging the terms and employing the redefinition of the coupling, A — 7r/|A|, 
one obtains 

^ = -'5^ - 2^0 iv, M, T) + ASv{v-<Jc) Ii {v,fi,T)-6^v' Pi iv, H,T)r - S^i^v, A^, T)]' , (4) 

where 



(T) AT) , AT) p 



Setting 6=1 and applying the PMS to Eq. ^ we obtain a self-consistent equation for the optimized mass 



o-c + ^riiXi{'n,ix,T)+l2{Tl,^.,T)I'^{'n,ii,T) [li{-q, ii,T) + ■qX[{ri, ^i,T)] ^| _ , (6) 



where the primes indicate derivatives with respect to rj. It can easily be verified that one immediate consequence of 
Eq. ([5]) is that for — > 00 we have fj = ac, thus retrieving all large- A^ results, as mentioned in the previous paragraph. 
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All the details concerning the Matsubara sums and dimensional regularized integrals can be found in Ref. fl^ . (Note 
that the sum over the Matsubara frequencies and all space momentum integrals can be evaluated analytically at this 
first S order here considered). 

In the case fi — 0, Eq. ^ factorizes in a nice way which allows us to understand how the OPT resums the series 
producing non-perturbative results. Using l2{r],0,T) = (see Ref. we can write Eq. ([6]) in the self-consistent 

form: 

fi = <J, + ^fjI,{fj,0,T) . (7) 
Explicitly, Ii{fj,iJ,,T) is given by (from the Matsubara sum and momentum integration) 
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The resummation becomes even easier to be visualized if one notes that Eoxc(^, 0,T) = Trfj/{NA) Ti{fj,0,T), 
where Eoxc(^7 0, T) represents an exchange (rainbow, Fock) type of contribution so that, at fi = 0, one may write 
fj ~ <Jc + Scxc(?7, 0, T). Note that prior to optimization, the terms contributing to VeS to first order are the ones given 
by Eq. ([3]), represented by the diagrams shown in Fig. [TJ Upon optimization the diagrams are dressed by rainbow 
type of self-energy contributions that originates from the first nontrivial diagram represented by the third term in Fig. 
[TJ This was expected since at order-J the perturbative OPT effective potential receives information about this type 
of topology only. The OPT will then gather all powers of corresponding to the (rainbow) class of diagram. In 
other words, the simple evaluation of a first topologically distinct diagram is already able to bring non-perturbative 
information concerning that type of contribution. If one proceeds to order-(5^, information about corrections to the 
scalar propagator as well as Yukawa vertex will enter the perturbative effective potential. Then, the PMS will dress 
up these perturbative contributions and so on. Of course this rather simple OPT/PMS procedure does not give all 
the truly non-perturbative and higher order contributions, such as the genuine and higher order corrections, but 

is just one particular class of graphs that are being resummed in this way. It nevertheless gives a non trivial result 
that goes beyond the standard 1/A'' order. 

Let us now consider the T = fi = case. By optimizing V^ff (uc), at 0{S), it is very easy to show that one finds the 
optimum value fj — acJ'{N), where CTc is the minimum of Vcs{a'c), and J-{N) = 1 — 1/ (4A^) is a characteristic function 
that will appear in all our expressions. At this value of 77 and computing the cxtremum for the effective potential for 
CTc, we obtain that the order parameter is given by 



<Jc = AT{N)-^ . (9) 

Note that the scalar field vacuum expectation value (Tc sets the natural scale for CSB. We could express all our other 
results in terms of CTc, but here we find more convenient to normalize all critical quantities in terms of just A (after 
all this is just a scaling choice). In Ref. the result for the critical temperature has been obtained in a similar and 
completely analytic fashion producing the OPT result 



T,{^l = 0)^A/[2\n2T{N)] . (10) 

However, an interesting alternative is to obtain Tc using the Landau expansion for the free energy density, which is 
valid for small values of the order parameter. The Landau's expansion (for small Uc) for Ves can be expressed in the 
general form 



Kff(fTc, T)^Vo + ^a(/i, T) al + h,{ii, T) at + \c{ti, T) al , (11) 

where Vq is a constant (field independent) energy term. Note that only even powers of are allowed due to the 
original chiral symmetry of the model. The coefficients a, h and c appearing in Eq. pi|) can be obtained, respectively, 
by a second, fourth and sixth derivative of the free energy density expansion at tJc = 0- Higher order terms in the 
expansion ([TT|) can be verified to be much smaller than the first order terms and can be consistently neglected. 

Lets recall that the expression of the free energy in the form (|11[) is the simplest one able to exhibit a rich phase 
structure, allowing also for studying tricritical phenomena. In particular note that a tricritical point can emerge 
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N 


Tc,Eq.^ 


rriPMSl 

C 


^PMS2 


1 


0.961797 


0.979106 


0.963579 


3 


0.786925 


0.786983 


0.786925 


4 


0.769438 


0.769454 


0.769437 


6 


0.752710 


0.752713 


0.752710 


10 


0.739844 


0.739844 


0.739844 



TABLE I: The critical temperature, Tc (at /i = and in units of A), obtained with the OPT in three alternative ways for 
different values of A'^. We have considered as much digits as we needed to establish eventual numerical differences. 

whenever we have three phases coexisting simultaneously. For the above potential, Eq. (|lip . one can have a second 
order phase transition when the coefficient of the quadratic term vanishes (a = 0) and 6 > 0, c > 0. A first order 
transition happens for the case of & < 0, c > 0. The tricritical point occurs when both the quadratic and quartic 
coefficients vanish, a = b = (with c > 0). Thus, Eq. (jlip offers a simple and immediate way for analyzing the phase 
structure of our model. For instance, to obtain at /i- = one only needs to consider Eq. (fTT|) to order-cr^ with 6 > 
to assure that the potential is bounded from below. Then, the solution of a(0, Tc) — sets the critical temperature. 
However, in order to use Landau's expansion we must have VcB in terms of ctc, /i and T only (apart from N and the 
scale A, of course). In principle, this can be done by using the PMS relation, Eq. ([6]). As one can easily check the 
large N result, Tc — A/(21n2) is quickly recovered by solving a{0,Tc) — with the latter obtained from Eq. ^ with 
fj = Uc- Now, at finite N , fj depends on CTc in a highly nonlinear way and the use of Landau's method in conjunction 
with the OPT-PMS does not appear to be straightforward. However, Eq. ([7]) can then be easily solved numerically 
by iteration. Let us consider the first order (in the iteration) approximate PMS solution (PMSl) 

(12) 

and to access the convergence of the solutions, let us also consider the next approximation (PMS2) by doing fj fj^-^'' 
on the right hand side of Eq. ([7]). In Table I we compare the results given by PMSl and PMS2 with the analytical 
relation Tc = A/[21n2 J^(iV)], Eq. (dU]), for /x = 0. 

This shows that even the PMSl approximation gives excellent results supporting the possibility of applying the 
OPT-PMS in conjunction with Landau's expansion. Especially, we notice that the results quickly converge for > 4 
but even for the extreme case A^ = 1 the agreement is rather good (higher than 99%). 

The reader should carefully note that we have written our results, in Table I, up to seven digits just to compare the 
two numerical iterative optimization procedures, PMSl and PMS2 with the analytical solution, Eq. (|10p . We have 
considered as much digits as were needed in order to establish the eventual numerical differences. By no means this 
otherwise excessive number of digits is intended to refiect the intrinsic accuracy of the OPT method at lowest order. 
At the end of the next section we briefly comment on the accuracy of the OPT based on the results at next order 

IV. THE PHASE DIAGRAM IN THE OPT 

Let us now present the phase diagram for the three dimensional model. Since the results for the cases T ~ ^ — 
and r 7^ 0, ^ = have been discussed in the previous section let us discuss the case T ~ and ^ ^ for which the 
relation fj = Uc ^{N) is still valid (the reader is referred to Ref. [3| for details concerning the evaluation of integrals 
and Matsubara's sums). In this case, CSR is found to happen at a critical value for the chemical potential, /ic, which 
is obtained as follows. We first note that for fj, = (and T = 0, as shown previously) the effective potential has 
a minimum at the value of ac = AJ^{N)~'^, as evaluated above. For 7^ the effective potential acquires another 
minimum that gets degenerate with the first at a value of chemical potential fic, where CTc = 0. At this point we then 
find that T^ff (ctc = 0, /i = /ic, T = 0) = Vcs{d'c, = 0, T = 0), from where we get, after some algebra, that the critical 
value for the chemical potential satisfies the equation 

Eq. (|13p is a forth order equation for /Xc/A, which is much simpler to solve by iteration, converging very quickly since 
fi/A is 0(1). For example, at A^ = 3, one finds Hc — 1.067A, which agrees with a full numerical evaluation of the 



fj^'^ ^^c + ^fjIiifj,0,T) 
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complete T and /i dependent expression ([3]) at 0{5) Since there is a potential barrier to overcome at the value 
of /i = the chiral transition here is discontinuous, thus first order. Note also that, for N — > oo, Eq. gives 
lie ^ and we recover again a well known large- result 0. 

Finally, we analyze the finite T and /i effective potential. This is done numerically for Eq. ^ at 0{6). For each 
value of the temperature and chemical potential we apply the PMS optimization condition to Vcs determining the 
optimized value of rj, fj = fj{(7c, This is inserted back in the expression for V^ff , from where we then obtain the 

regions for CSB and CSR. The result for the phase diagram, for = 3, is shown in Fig. [21 




FIG. 2: The phase diagram for N — 3 obtained within the OPT. The position of the tricritical point is indicated by Ptcr, where 
the second order transition line (above Ptcr) turns into a first order one (below Ptcr). CSB happens at the left of the continuous 
curve while CSR happens at the right hand side. The dashed and dot-dashed lines are the metastable lines running along the 
first order line. The inner plot zooms into the region around the first order transition line. 



According to the large-A^ approximation the discrete chiral transition in the GN model is of the second order 
everywhere except at T = and fi = fi^- The results obtained from the OPT, however, show the existence of a 
tricritical point at Ttcr — 0.25A and ^tcr — 1.03A (for TV = 3), represented by the black dot in Fig. [21 Below this point 
the transition is of the first kind, while above it is of the second kind. Metastable lines can be found at the left and 
right of the first order critical curve on Fig. [21 running very close to it from the tricritical point down to the fj, axis. 
Any method relying on determining this region, if not enough precise, is very likely to miss it. This may then possibly 
explain why previous works, most notably lattice Monte-Carlo simulations 5], had difhculties to probe the possible 
existence of a mixed phase. It is possible to use the OPT-PMS in conjunction with Landau's expansion in order to 
determine the tricritical points for any finite values of N. In this case one can just reinstate the chemical potential 
back in the last term of Eq. (O, iterate it to the desired order, substitute the solution back in the expression for the 
effective potential to then obtain the coefficients a(/z, T) and 6(/x, T) in Eq. (fTTjl . Then by solving a{fj,tcr, Ttcr) = and 
b{fitcr,Ttcr) = simultaneously, the tricritical points can easily be extracted. As discussed above, direct numerical 
evaluation using the Eq. ^ and the PMS optimization give Ttcr — 0.25A and /iter — 1.03A for = 3. The Landau's 
expansion, with for example the 0PT-PMS2, gives Ttcr — 0.24A and /iter — f .03A. These results again compare very 
well with the full numerical results used to obtain Fig. [21 In Table II we present all the relevant results for chiral 
symmetry restoration for different values of N. The quantities Tc (at /i = 0), Ttcr and fitcr have been originally 
obtained in the present work using the Landau expansion within the next order iteration PMS2, as described above, 
while fic (at T = 0) is given by the solution of Eq. (|I3p and the order parameter, CTc, is given by Eq. 

It is also possible to express and analyze the phase diagram in terms of the pressure P and density p, defined 
respectively by P = —Vcs{v,o'c,l^,T) and p — dP/dp,. This was done in [l3|, where it was shown the existence of 
a chiral broken/restored coexistence phase (the analog of a liquid-gas phase), whose result is a direct consequence 
of the existence of the tricritical point seen in Fig. [21 At this point one may wonder how the consideration of the 
order-5^ terms, which contains additional \/N as well as new corrections, would affect the first order results. 

Clearly, at T 7^ and p the explicit evaluation of all the three loop contributions is a rather difficult task. 
Nevertheless, at T = and p = 0, a. lengthy evaluation of those contributions has been done in Ref. showing 
that the potential corrections to the 0{S) results are not negligible but remain under control, about a few percent. 
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0) 


Ttcr 


/iter 


1 


1.78 


0.96 


1.24 




0.46 


1.06 


3 


1.19 


0.79 


1.07 




0.24 


1.034 


4 


1.14 


0.77 


1.05 




0.21 


1.03 


6 


1.09 


0.75 


1.03 




0.19 


1.02 


10 


1.05 


0.74 


1.02 




0.17 


1.01 


oo 


1.0 


0.72 


1.0 




0.0 


1.0 



TABLE II: The relevant physical quantities (in units of A) for different values of A*'. 

For example, the result obtained for the scalar field vacuum expectation value at T = and /i = is found to be (at 
TV = 3) (Tc ~ 1.17A at 0(5^), for a reasonable choice of the arbitrary renormalization scale, chosen as given by A, 
while at 0{5) it is CTc — 1.19 A. For larger values of N the difference gets even smaller, in conformity with the large- A'' 
results (since at A^ — > oo the OPT converges exactly to the large- A'^ results [IH). This then hints on the accuracy and 
quick convergence of the order-5 results discussed in this letter. 

V. CONCLUSIONS 

In conclusion, by studying the 2-1-1 dimensional four-fermion GN model with discrete chiral symmetry in terms of 
the OPT, we have shown how this method can be advantageous compared to the widely used large- A^ and Monte Carlo 
approximations. Here we have presented precise analytical expressions for the order parameter CTc (at T = fj, = 0), 
the critical temperature Tc (at /i = 0), and the critical chemical potential fic (at T = 0) for GSR which are valid 
for any values of A^ already at lowest order in the OPT. As usual, all these relations exactly reproduce the standard 
large- A^ results when the N ^ oo limit is taken. We were also able to precisely locate a tricritical point in the (T, /i) 
phase diagram not seen using the large- A^ approximation and also not obtained before using other methods, though 
numerical simulations seemed to indicate its presence. We have discovered that, due to the finite A^ effects, the model 
behaves as a Van der Waals liquid in which a coexistent liquid-gas phase does exist. The observed metastable region is 
rather small, which possibly explains why its observation was missed in Ref. [EJ and previous works. All these results 
seem to survive when going to higher orders in the OPT, as we discussed in the last section. In this letter we have 
extended our previous applications of the OPT to the Gross-Neveu model [U [3] by analyzing how the method can 
be applied to the Landau expansion for the free energy density in order to determine the value of critical quantities 
in a way that is less time consuming than the straight numerical methods used in Refs. [HI, [l3|. We have shown that 
this is possible if one considers further approximations in the optimization equations. However, the numerical effects 
of such approximations, even at the first level, turn out to be negligible. 

Goncerning the next order OPT corrections we have checked in Ref. [13] that, for T = = ii, N — 3, and a 
reasonable choice of renormalization scale the differences for the order parameter are about 2%. As for the 1/A^ 
expansion we are not in position to speculate what a complete next to the leading 1/A^^ order calculation would give 
from a quantitatively perspective. However, we believe that it should point out to the existence of a tricritical point in 
the T — fi plane since, according to our work, this seems to be an effect of finite A^ corrections which is also supported 
by Monte Garlo applications [JQ. 

Possible follow-ups to this work that would be interesting to investigate regard the possible existence of a kind 
of nuclear matter state around fic-, at large fi and low T, for which the chiral phase would be degenerate with a 
baryonic phase, (ipip) and {ip^otp) 0- This is still an unsettled question in the literature and that could be 
ideally addressed with the OPT. Another interesting work would be to evaluate the various critical (and tricritical) 
exponents, determining whether they belong to the mean-field universality class or the Ising universality class in two 
dimensions. This is also a controversial point concerning the GN model, with no clear answer so far in the literature. 
Based on its recent success in dealing with critical theories [ll|, HI, Hli we expect that the method employed here will 
be a valuable tool in these analysis providing us with a better understanding of the chiral symmetry behavior in a 
hot and dense medium. 

Finally, we point out that since the GN model is part of a larger class of four-Fermi models, e.g. the Nambu- 
Jona-Lasinio type of models commonly used in nuclear physics problems and as an effective model for QCD, or even 
its non-relativistic version, used in the study of planar superconductors and low dimensional fermionic systems (like 
Fermi gases in atomic physics), we expect that not only our results but also the OPT method used to obtain them 
can be of interest and have applications in these areas as well. In particular, we believe that the success of the OPT 
method in redrawing the phase diagram of the three dimensional GN model could certainly attract the attention of 



8 



condensed matter as well as nuclear physicists to this method as an alternative tool to the well known large- iV and 
Monte Carlo techniques, which are currently used in an interdisciplinary manner. 
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